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Abstract 

The behavior of a thin him of nematic liquid crystal with unequal Frank con- 
stants is discussed. Distinct Frank constants are found to imply unequal core 
energies for +1/2 and —1/2 disclinations. Even so, a topological constraint 
is shown to ensure that the bulk densities of the two types of disclinations 
are the same. For a system with free boundary conditions, such as a liquid 
membrane, unequal core energies simply renormalize the Gaussian rigidity 
and line tension. 
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I. INTRODUCTION 



This paper discusses the disclination-mediated isotropic-ordered transition in a thin film 
of nematic liquid crystal. The focus is on the case where the bend and splay Frank constants 
are distinct. The free energy of this system is given to lowest order by [I 

H = y J dxdy [V • n(x, y)f + y J dxdy | V x n(x, y) | 2 . (1) 

Here n is the orientation of the nematic molecule and is of unit length. This free energy can 
alternatively be expressed in terms of the orientation of the molecules as 

H = ^J dxdy (e 2 x ± 9f) ± — J dxdy cos(2#) (flj - Of) 

-A j dxdysm(26)8 x 6 y , (2) 

where J = (k\ + ks)/2, A = (hi — ks)/2, and the subscripts denote derivatives. If the two 
Frank constants were equal, this free energy would simply be that of the X-Y model. For the 
symmetric nematics considered here, however, the natural defects are ±1/2 disclinations, 
rather than the ±1 disclinations of the conventional X-Y model. It will be shown that the 
presence of non-zero A causes the ground state energies of ±1/2 disclinations to differ by 
0(A 2 ). The disclination energies diverge logarithmically in the system size, but with unequal 
coefficients. The elementary Kosterlitz-Thouless energy-entropy balance thus seems to lead 
to different proliferation temperatures for these defects. This famous argument @ predicts 
that a ±1/2 or —1/2 defect proliferates whenever the free energy to create a disclination, 
F +1/2 (R) = E +1/2 (R) - 2k B T\n(R/a ) or F_ 1/2 (i?) = E- 1/2 (R) - 2k B T\n(R/a ), becomes 
negative. Here E + \/2{R) and E-\/2{R) are disclination energies as a function of the system 
size R, and ao is a microscopic cutoff. 

In fact, thermal fluctuations of the nematics drive the two Frank constants to the same 
value at long wavelengths, so that there is a unique Kosterlitz-Thouless transition temper- 
ature. The essential effect of non-zero A is to create a distinct long-ranged contribution 
to the core energy of each defect. For a system above the isotropic-ordered transition, this 
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difference in core energies can be substantial and would seem to lead to different densities 
of +1/2 and —1/2 disclinations. As the correlation length grows near the isotropic-ordered 
transition, the core energy becomes negligible compared to the logarithmically diverging 
piece, and the disclinations pair into dislocations. At the transition, the densities of +1/2 
and —1/2 disclinations must both become equal to the density of dislocations. In fact, 
Green's theorem implies that under all conditions the difference between the number of 
+1/2 and —1/2 disclinations can scale at most as the circumference of the system. The 
natural way for this to happen is n+1/2 — W-1/2 ~ c -R/£> where £ is the correlation length, 
and the prefactor depends on difference of exponentials of core energies. This constraint 
implies that ±1/2 disclinations occur with the same density in a large system. 

The ground state energies of ±1/2 disclinations are derived in Sec. II. The energies 
are found to be logarithmically diverging with unequal prefactors. It is shown in Sec. Ill 
that these coefficients should renormalize to the same value at long wavelengths due to 
renormalization of A to zero. The free energies of ±1/2 disclinations in the ordered phase 
are directly calculated by perturbation theory in Sec. IV, and the ±1/2 disclinations are 
found to differ by a core energy contribution. An approximate calculation of the disclination 
density in the isotropic phase is described in Sec. V. The ±1/2 disclinations are shown to 
occur with equal densities for large systems. In fact, the number asymmetry is shown to 
scale only linearly with the system size. This issue is explored in Sec. VI with Monte Carlo 
calculations on a lattice model. The disclination number asymmetry is indeed found to be 
proportional to the circumference of the system. Section VII concludes with a discussion of 
these results. 

II. GROUND STATE ENERGIES OF ±1/2 DISCLINATIONS 

The ground state energies of ±1/2 disclinations in the Hamiltonian (|2]) will be loga- 
rithmically diverging in the size of the system. The coefficient in front of the logarithm 
is calculated in this section. The coefficient is determined solely by the properties of the 
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9 field far from the disclination, vanishing if the 9 field vanishes at infinity and diverging 
unless the 9 field goes to a constant. The ground state configuration, therefore, must have 
9(r, 0) ~ 9 ((f>) as r — > oo. With this form, the energy is given by 

H ~ i7 M-R/oo) as -> oo , (3) 

with 

H = - d<p 9' 2 (<P) { J + A cos[20 - 26(<f>)]} . (4) 
2 Jo 

A disclination of strength s is defined by 

8(<f>) = s<j>+j8 1 (<i>) , (5) 

where 9\ (0) is continuous. The condition of a ground state geometry, without the assumption 
of rotational symmetry, can be written as 

y = = -JV 2 £ + Asin(2(9)(e; - 9\ + 29 xy ) 

+Acos{29){9 xx - 9 yy + 29 x 9 y ) . (6) 

This equation implies that 

9 = 9 + j9 1 + O(A 2 /J 2 ) , (7) 

with 

V 2 9 = 

V 2 #! = sin(2^ )(^ " C + 2^,) 

+ cos(2#o)(#cte:r — 9o yy + 29o x 9oy) . (8) 



These equations are solved by 
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for a disclination of strength s. 

For s = +1/2, Eq. (0) simplifies to 



and leads to an energy of 



0i(0) = -|bM (10) 



Similarly, for s = —1/2, Eq. (H) simplifies to 



and leads to an energy of 



^(0) = ^sin(30) (12) 

60 



The general result for the deviation field 6 — 6 can be expressed in terms of elliptic 
integrals of the third kind ||. Results for the angle field have been presented elsewhere [|J. 
To give the reader some feel for how larger values of A distort the ground state geometry, 
the function 9(4>) is shown for a s = +1/2 (Fig. 1) and a s = —1/2 (Fig. 2) disclination. 
These geometries were computed by defining 6 ((f)) on a grid and numerically minimizing Eq. 
(f|). Direct integration of Eq. (^) produced identical results. Extreme differences between 
the two Frank constants can substantially distort the T = geometries. Figure 3 shows the 
ground state energies associated with different ratios of the Frank constants for s = ±1/2. 
One can see that the s = +1/2 disclination completely screens out either splay or bend as 
the associated Frank constant, k\ or k 3 , respectively, becomes large: 

H ~ as ki — > oo 

H ~ nki/2 as k 3 — > oo (14) 

The s — —1/2 disclination, on the other hand, is unable to completely remove unfavorable 
bend or splay: 

H ~ 0.191&1 as k\ — > oo 

H ~ 0.191Jfe 3 as k 3 -> oo (15) 
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III. RG FLOW EQUATION FOR A 



The isotropic-ordered transition does not occur at T = 0, and so it is the free energies of 
the two disclinations that should govern their densities in the isotropic phase. This issue is 
addressed here by looking at the renormalization of A due to thermal fluctuations around 
an isolated disclination. 

For a s = +1/2 disclination the angle order parameter is expressed as 

9(r, 4>) = <f>/2 + V»(r, <f>) . (16) 

The function if)(r,<j)) is single- valued and smooth. The Hamiltonian ($) is expanded in 
powers of ip, with the result 

l- _ . n i .-1 



a i- 
+ lJ 


3x?p 2 | 2x(^J- 
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3y^ 3 Qy^y- 
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A f 
+ T2l 






r 



(17) 

The renormalization of A is tracked by integrating out i/j on a momentum shell to first 
order in T/J and watching how terms from the ip 3 expression contribute to terms in the ip 
expression. The ip expression can be written as 

= — lf* ( - k) - (18) 

where / k means / d 2 k/(2ir) 2 . The if} 3 expression can be broken down into 

t! = -TriA[ (27r) 2 (5(k 1 +k 2 + k3 + k 4 )^ 

Jkik 2 k 3 k4 k\ 

X ^(k 2 )^(k 3 )^(k 4 ) 

t 2 = —2mA [ (2n) 2 5(k 1 + k 2 + k 3 + k 4 ) 

ikik2k3k4 k\ 
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x ^(k 2 )^(k3)^(k4)(/c 3a; fc4x - hyk. 



■hi) 



-47iiA 



(27T) 2 5(k 

kik2k3k4 

x ^(k 2 )^(k 3 )^(k 4 )(-/i;3 ;c A;4j / ) . 



k 2 + k 3 + k 4 )-^ 



(19) 



The momenta in the shell k c — dk c < k < k c can then be integrated over with the result 

u \ -a f k v 1/ i x67rTfc c dfc c 



(*2>, 








This result implies 



A' = A 



AT rffc. 



7rJ k c 

Defining dk c /k c = dl, the flow equation results: 

dA AT 



dl n J 

J is not renormalized to 0(A). 

The same calculation can be performed for a s = 

0{r,4>) = -4>/2 + i[)(r, 



(20) 



(21) 



(22) 



T/2 disclination, using the relation 



(23) 



The Hamiltonian @ is expanded in powers of ip, with the result 
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5(a; 3 — 3xy 2 )ip 2 


+ Jj 
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+ ll 


5(3x 2 y — y 3 )ip 3 




A /■ 
+ I2/ 


5(3xy 2 — x 3 )ip 4 
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y • x' _|_ 
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(24) 

The renormalization of A is again tracked by integrating out t^ona momentum shell to first 
order in T/J and watching how terms from the ip 3 expression contribute to terms in the ip 
expression. The flow equation that results is 



dA AT 

n = -7f ■ < 25 » 

One can see that A renormalizes in the same way about ±1/2 disclinations. In fact, the 
same flow equation describes the renormalization of A in the absence of disclinations ||. 

The renormalization of coupling, J, and the disclination fugacity, y, can be studied with 
the correlation function approach used for the standard X-Y model 0. The result, combined 
with the results above, is 

^ = -AT/(nJ)+0(A 2 ,y 2 A) 

^=(2-^)^ + 0(, 2 A,, 3 ). (26) 

One can see that at the critical point the renormalized coupling is Jr/T = 8/tt. Furthermore, 
A scales as 

AR-A^/ao)" 1 / 4 , (27) 
where the correlation length is given by £ = a exp(Z). 

IV. DISCLINATION FREE ENERGIES IN ORDERED PHASE 

While the RG calculations show how A becomes irrelevant at the isotropic-ordered tran- 
sition, they do not directly show how the disclination densities scale near the transition. To 
estimate the ±1/2 disclination densities, expressions are needed for the free energies of iso- 
lated disclinations. Perturbation theory is here used to calculate directly these free energies 
in the ordered phase. The results should also be applicable to correlated regions within a 
macroscopically disordered phase near the isotropic-ordered transition. 

The free energy of a disclination will be evaluated to 0(A 2 ). Equation @ is first inte- 
grated by parts with the result 

J f . . /„o „o\ A 



H=-f dxdy {el + 9 2 y )~-J dxdysm(26) (9 yy - 6 XX ) 

[ dxdycos(29)6 xy . (28) 



Although 9 is discontinuous in the presence of a disclination, the result is the same if Eq. 
(|T6|) or (|23| ) is used and the integration by parts done in terms of if). A cumulant expansion 
is used for the free energy: 

F = E + (6H) 0c - ~ ((^) 2 ) 0c /T + ... . (29) 

Here the averages are done with respect to the reference system with Hamiltonian H Q , 
indicated by subscript zero, and are connected, indicated by subscript c. The functional 5H 
is H — H . The reference system is chosen to be 

W] = |/dr(^a+^ + ^) , (30) 

where 9 = scf) + if>. 

Specializing to the case of a s = +1/2 disclination, the perturbation becomes 



5H\ib\ = / dr 



— - sin [2w) 



sin(2^) - — cos(2ip) (if> yy - ip xx ) 



-— cos(2if))if; xy + cos(2^) 

-— sin(2ip)(i/)yy - ip xx ) + — &n(2ip)ip xy 



(31) 



A short calculation shows 



(5H) 0c = . (32) 

The first non-zero contribution to the free energy is, therefore, 0(A 2 ). The form of Eq. 

with three terms even in if) and three terms odd, simplifies the evaluation of ((5H) 2 ) 0c . 
Even so, there are twelve Gaussian averages that must be performed. A typical term is 

J dr 1 dr 2 (cos[2^(ri)]^ xa; (ri) cos[2^(r 2 )]Vw( r 2)) 0c 
= 2 J dr lC Zr 2 e 4 * (ri2) - 4 * (0) 

x{[Xxx{rn) - Xxx{0)][Xxy{ri2) ~ Xxy{0)} + ^Xxxxy(r 12 ) 
4%(ri 2 )-4x(0) 



+ 2 J dridr 2 e 



X\ - \Xxx( r w) + Xxx(fy}[Xxy(ri 2 ) + Xxy(0)} + ^Xxxxy(r 12 ) 



(33) 



where T\ 2 — Irj — r2|, and 



T 

m 



jk 2 

T 

X{r) = ln[R/r) 
In J 

X (0) = -^-ln(«/oo) • (34) 

With the definition z = 2T/(ttJ), the first integral in Eq. ( P5| ) scales like (R/a )~ z , and 
the second term scales like (R/ao)~ 2z (R/ao) z Both terms must in principle be evaluated. 
However all twelve terms that contain the factor e~ 4x ^' 12 ^~ ix ^ cancel by x «-» y symmetry. 
The symmetry r x <-> r 2 is applied to the other twelve terms with the result 



| dr 1 dr 2 ((6H) 2 ) Qc = 

A 2 r 
+ — / dr 1 dr 2 e Ax( - ri2) 
16 J 



+ <{ 4— — [x ra (r 12 ) - X^(ri2)] 2 + 16— —xl y {ni) 



+ {^V*x(r 12 )) 
In r 2 J 



.ri r 2 

+ {-4^^[x„M-^(r 12 )]-8^^ (ris )' 
[ r l r 2 r l r 2 

= h + 7 2 + 7 3 + h (35) 

with the redefinition xij'12) = ~ (-2/4) ln(r 12 /a ). The four integrals, I\ — J 4 , represent inte- 
gration over the four terms in braces. 

These integrals can now be evaluated. This will be done in Fourier space, and the 
following Fourier transforms will be helpful 

I J 2 z T(z/2) 
FT^y/r 3 } = 2m sin 

FT{y/r} =2m sin <f)/k 2 . (36) 
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The first of these transforms is well-defined for < z < 2. The relation is valid for all z by 
analytic continuation from the relation FT{V 2 f} = —k 2 f(k). With these definitions, the 
first integral becomes 

A 2 a* 



h = 



16 

l 2^2 



/ |2«sin0| 2 2vrF- 2 ^-4^ 
k 1 Yl 2 z T{z/2) 



A 2 tt 2 2r(l - z/2) 1 - (R/a )- z 
~~ 8 2 z Y(z~/2) z 
A 2 vr 2 2 T(l - z/2)(l + z/2) 1 - (R/a y 



8 2^r(2 + z/2) 



The second integral becomes 



lb J \a / r\ r 2 

2 



16 



2iri sin 



k 2 



2rf+2 2T(-1 - ,/2) 
2 4+2 r(2 + ^/2) 
_ A 2 ttV 2r(-l - z/2) 1 - (i?/a )- 2 
8 2 4 +*r(2 + 2/2) z 
A 2 ti 2 z T(l-z/2) l-{R/a )- z 
~ 32 2 Z T(2 + z/2)(l + z/2) z ' 

The following identity will be useful in evaluating the third integral: 

/ \ z 1 - (r/a o y a 

X(r) = -- lim , 

4 Q -*o a 



(37) 



(38) 



(39) 



which implies 



V\(r) = -lima(a + 2) 2 (r/a )- a - 4 



(40) 



Using this result, one finds 



z- l FT{r- z V\\ =-\ima(a + 2) 2 - — — 

I A I 4«^o v > ( z + a ) 2 (z + a + 2) 2 



x 



= < 



47ik z+a - 2 T[l -(z + a)/2] 
2 z + a T[(z + a)/2] 

nk 2 /2, z = 
0,2 ^ 



(41) 
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The term J 3 , therefore, vanishes for non-zero temperatures. It is convenient to break the 
fourth integral into two parts: 



ha = I dridr 2 ^— (—) [Xxx{r 12 ) - X w (ri 2 )] 
4 J r±r 2 \ao / 

^46 = — — / dritfo-^ — ( ) Xxi/(ri 2 ) • (42) 



>? r 2 V a 

— / arictr 2 ^- — — 

2 J rfr 2 \a 

The following trick is used to evaluate these terms: 



r- z [Xxx(0 - x ro (0] = Mm -r~ 2 (^ - d 2 y )r~ a 



4a 

V*- 2 cos(20) . (43) 



In this form, the Fourier transforms can be evaluated, with the result 



ha — 



,2^2 



^/^(-?)^^i 



Similarly, 



2^r(2 + z/2) 

A 2 ti 2 z T(l-z/2) l-{R/a )- z 
16 2 2 r(2 + ^/2) z 



A 2 Oq r . 27ricos0 / nz 



hb = ^Jj2mam<j>) 

x sin(20)/c 



A; 2 V 4 

m-z/2) 



2 z T(2 + z/2) 
A 2 n 2 z F(l-z/2) l-(i2/oo)- 



(44) 



(45) 



16 2 2 r(2 + ^/2) z 
Combining all these results one finds for the free energy of a +1/2 disclination at the origin 

F +1/2 = ^-ln(R/a ) 

ttA 2 Hi -z/2) / n 1 \ 

v 1 ' 1 c + 2z + ' 



32J 2 2 r(2 + ^/2) V 1 + z/2 J 

x l-W°o)-' +0{A s h (46) 



The ground state energy is found to be 

/ vrJ 9ttA 2N 



F +1/2 ~ ^ — __ j \ n (R/a ) + 0(A J ) as T - , (47) 
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in agreement with Eq. (|TT|). 

For the case of a s = - 1/2 disclination, the perturbation to consider is 



A 



5H[ip] = — I r/r 



r y 3 — 3x 2 y . 



y 



sin(2^) + - cos(2ip) (ipy y - ifj xx ) 
cos(2ijj)ip X y H cos(2?/>) 



r 
x 



.1 



— sm(2ip)ip X y 



(48) 



A short calculation shows 



<«Ooc = , 



(49) 



so that the first non-zero contribution to the free energy is 0(A 2 ). As before, averages 
that lead to terms with e~ 4x(ri2 ) _4x< - ' ) cancel by x <-» y symmetry. Also as before, the term 
containing the V 4 x{ r ) vanishes at non-zero temperature. After some simplification, one 
finds 

A 2 



Oc 



where fi = (yf — 3x\yi)jr\. The result 



/2 + T6 J dr ^ e4x(ri2) [ifM 

+ {±h — [Xyy{n2) ~ Xxx(r 12 )} 

+ \-8fi—Xxy{ri2, 
I r 2 
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(50) 



/(k) = vk ( —6 sin + 8 cos 2 (j) sin <\> + — sin 3 1 



will be used. 

The integral ([50]) is split into the three bracketed pieces. The first integral is 



(51) 



AV 



1 1 |/(k)| 2 27TF- 22r(1 "- /2) 



16 A"""" 2 2 r(^/2) 
A 2 7r 2 ^r(l - z/2)(l + z/2) 1 - (R/a y 



72 



2*r(2 + z/2) 



(52) 



The second integral is 
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— / dridr 2 — 4/i — [x w (ri 2 ) - Xxx^u)] 
16 j V ao / ?"2 

A^ f /t(k) 2^in^ cos(2 ^ r(l-,/2) 
4 «/ k A/ 2 



4 ik w k 2 2 v ^ 2 z r(2 + ^/2) 
A 2 tt 2 ^ T{l-z/2) l-{R/a )- z 
48 2 2 r(2 + z/2) z ' 



(53) 



The third integral is 



J 7 = -/«r 2 (^) Z ^8 fl ^) Xxy ( 
A 2 ag f ~ 27ricos0 / ttz 



r 12, 



2 ik J w k 2 V 4 

x s in(20)F r( !~" :/2) , 
v V! 2 z T(2 + z/2) 

A 2 n 2 z T(l-2/2) l-(R/aoY 



(54) 



48 2 z r(2 + z/2) 

Combining all these results one finds for the free energy of a —1/2 disclination at the origin 

F_ 1/2 = ^-\n(R/a ) 

vrA 2 r(l-z/2) / 9 \ 

'16 + 2,2 + 1 



288J 2^r(2 + z/2) \ 1 + z/2 J 

x l-(iVa )-* +o(A3) - (55) 

The ground state energy is given by 

F_ 1/2 ~ ^ - ?^£pj HR/a Q ) + 0(A 3 ) as T - , (56) 

in agreement with Eq. fljjl). 

The difference in free energies of ±1/2 disclinations is, therefore, given by 

vrA 2 T(l-z/2) 7 + 2zl- (R/a )- z 



F+i/2 — F-1/2 



J 2 z T(2 + z/2) 36 z 
+ 0(A 3 ) . (57) 



Near the isotropic-ordered transition, the coupling renormalizes to zr = 1/4 and Jr/T 
8/71", so one finds 
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(F +1/2 - F +1/2 )/T ~ -0.8891(A/T) 2 [1 - ^/a )-^] 

+[E C+1/2 (A) - E C _ 1/2 (A)]/T 
as £ — > oo . 



(58) 



Additional microscopic core energies that may be distinct for the two different disclinations 
have been explicitly added in this equation. By comparing Eq. fl58|) with Eqs. flllD, flT3|) , 
and (p7|)) one sees that the appropriate finite-size scaling replacement is A^ln(i2/ao) — > 
A 2 [l — (i?/ao)~ z ]/-2- This relation gives a very good approximation to Eq. (|58|). 



V. DISCLINATION DENSITIES IN ISOTROPIC PHASE 

The disclination density in the isotropic phase is calculated in this section, taking into 
account interactions between disclinations. In the isotropic phase, the free energy can be 
expressed in terms of the order parameter 



Q = Qo 



Qo 
( 



^ ' n 2 x -n 2 y 2n x n y ^ 



V 2n x n y n 2 y -n 2 x ] 
^cos(2fl) sin(2#) ^ 



^ sin(2#) - cos(20) ) 
\ 



\Q2 -qi J 

An expression that reduces to Eq. (0) in the low-temperature, fixed-Qo phase is 



(59) 



ij ijk 

+ 32 J drQsiQtjdsdtQij , 



(60) 



where j\ = J/(AQq) and j 2 = A/(16Qo). Expressed in terms of the unique components q± 
and q2, the Hamiltonian becomes 



15 



H=jf dv(ql + q\) + || dv(q 2 lx + q\ y + q\ x + q\ y ) + 5H , (61) 

where 

5H = j 2 J dr q\(q lxx - q lyy ) + q\{q\ yy - q\ xx ) 
+2q 1 q 2 (q 2xx - <?2 ra ) 

+4qiq2qi xy + 2q\q 2xy - 2q\q 2xy . (62) 

The vector field q = (91,92) wm have disclinations of strength ±1 when the 9 field has 
disclinations of strength ±1/2. The density of disclinations can be written as 

P( r ) = ^5] s S n [ det, %(r)]<Hr - r,) , (63) 
where q(r^) = 0. This expression can be simplified as 

P( r ) = ^(l) det^-(r) 

= ^(9i)^(92)(9ix92^ - 9i^92x) • (64) 

Furthermore, the number density is given by 

|p( r )l = 2%0 5 (92) \qixQ2y - qi y q2x\ ■ (65) 

These densities will be evaluated by perturbation theory in A. The reference system 
will be Eq. flST|), with 5H = 0. It is clear that qi and 92 are independent Gaussian fields 
with correlation function Xo( r ) = T / (2irji)Ko(r j '£) , where the correlation length is given by 
£ = (ji/m) 1 / 2 . Averages that involve the factor S (91) will occur. This factor can essentially 
be absorbed into the Gaussian weight with the following trick: 

(5(a(0))f\a]) 

1 191 JJ;[9J ' JV[q]P[q} 

_ fV[q]P[q]6(q(0))f[q]fV[q]P[q]6(q(0)) 
JV[q]P[q}6(q(0)) J V[q}P[q] ' 

(66) 

where P[9]P[9] is the probability of a given field configuration. This result is defined to be 
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(/M) 5 (%(0))) = [27r Xo (0)]- 1 / 2 (/[g]) 5 . (67) 

It turns out that the weight P[q]S(q(0)) still implies that q is a Gaussian field, but with the 
new correlation function [pHTTIl 



(g(ri)g(r 2 ))«5 = x{ r i, r 2) = Xo(ki - r 2 |) 

- Xo(ri)xo(r 2 )/xo(0) . (68) 

The average number density of disclinations in the limit A — > is first evaluated. The 
following term will arise 

(5( qi (0))5(q 2 (0))\q lx (0)q 2y (0) - 9i»(0)g to (0)|) 

= pvrxoCO)]- 1 (\q lx (0)Q2v(0) ~ Qiy(0)Q2 X (0)\)s ■ (69) 

Using Eq. ( |d8|) , one can see that all of the variables in this average are independently 
Gaussian, with variance 7 = — Xoxx(fy = — Xoyyity = ~ Xor(O). The probability distribution 
for P(qi x q2y = x) is Kq(x / 7) / (^7) . And the probability distribution for P{q\ x q2 y — qi y q2x = 
y) is exp(— \y\ / j) / (2j) . With this result the average can be carried out: 

»iSf +0(A) - (70) 

The effect of non-zero A on the average disclination asymmetry is now evaluated. From 
Sec. IV, one knows that the essential effect of non-zero A is to create distinct core energies 
for +1/2 and —1/2 disclinations. This effect is modeled by replacing Eq. (|62j) with 

5H = -2 J dTfi(r)p(r) . (71) 

The parameter /x is related to the core energy difference, and there will be a unique mapping 
from A 2 to /i for small A. The core energy difference is considered to be non-zero and 
constant within a large region of radius R inside a macroscopically large system. This 
is done because the macroscopic disclination asymmetry is identically zero for a system 
with the periodic boundary conditions implied by the Fourier analysis. A grand canonical 
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ensemble with a fluctuating disclination asymmetry arises when considering a region of the 
periodic system. The disclination asymmetry is now calculated for small /x: 



(p(0)> = (p(0)) + 2im[ dr(p(0)p(r)) 0c + Ofa 

Jr<R 



(72) 



The disclination asymmetry is seen to be related to the disclination correlation function 
by a fluctuation-dissipation theorem. The disclination correlation function is given by || 



4<p(0)p(r)) = 2<|p|>o + <?o(r) 



(73) 



where the delta function comes from the self terms in the average, and go(r) is a radially 
symmetric function that accounts for correlations between distinct disclinations. To calculate 
go(r) the same trick as before is used. One has 

1 



w,(OM,(r)))o = Md(o) - 

The field correlation function is now given by ||TT 



(74) 



x{ r i, r 2) = Xo(|ri - r 2 |) - Xo{ri)xo(r)xo{r 2 ) 
-Xo(^i)Xo(^)Xo(|r 2 - r|) 
-Xo(\ri ~ T \)Xo(r)xo(rz) 
+Xo(|ri-r|)xo(r)xo(|r 2 -r|) 



(75) 



This form implies 



(?ia(0)gi/j(r))a = -Xoap(r) 

-Xo*(r)xo(r)xop(r)/ 1x1(0) - xi(r)} ■ 



(76) 



With these results one finds 



9o{r) 



^r[x 2 (0) ~ Xl(r)} 



X0r(r)X0rr(r) 



Xo(r)xlr( r ) 
Xm-Xl(r) 



(77) 
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Per Eq. ([Fj|) this result will be integrated over r. The first term can be integrated by parts, 
with the result 

This result, along with Eqs. ([72]) and Eq. (|75|), implies that the densities of +1/2 and —1/2 
disclinations remain equal even in the presence of distinct core energies. A more careful 
conclusion is that the average disclination density is zero within the bulk region. It can be 
non-zero near the boundary because the integral in Eq. (|72|) will be cut off before go (r) is 
negligible. This contribution is estimated to be 

(n+i/2 - n-u 2 ) = 2/2 drdr'(p(r)p(r')} + 0(n 2 ) 

Jr,r'<R 

« 27r//i?£<M>/2 

~ c[exp(-pE c+1/2 ) - exp(-/3£ c _ 1/2 )] 

xR/£ as £ < R -> oo . (79) 

The scaling (|p|) ~ c^ -2 has been used. It arises arises because the disorder created by 
unbound disclinations defines the correlation length. 

VI. MONTE CARLO CALCULATIONS ON A LATTICE MODEL 

This section describes both a lattice model for nematics with unequal Frank constants 
and a Monte Carlo procedure for evaluating the properties of the model. The fundamental 
degrees of freedom are spins of unit length on a square lattice. The Hamiltonian has both 
nearest- and next-nearest-neighbor couplings: 

i j<5 



1 - (rij ■ rij) 2 



2 | 2 _ 2 _ 2 

n ix ' U jx n iy U jy 



+fEE(-i) J 

° i i<5 
° i 4<j<9 
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(80) 



The sum over i is over sites on a square lattice. The sum over j is defined according to 
Fig. 4. Terms that would place j off the lattice are ignored. This Hamiltonian satisfies the 
symmetry rij <-> — iij. The couplings between spins are symmetric under the interchange 
rij <-> rij. Finally, in the limit of a very small lattice spacing, Eq. ( PPj ) reduces to Eq. (0). 
Numerical values of Eq. flgU| ) agree with those from Eqs. (0)-(§) for specific forms of the 
field. 

This model will be equilibrated with a simple Metropolis move that perturbs individual 
spins. Specifically, a vector randomly distributed in a disk of radius r is added to a randomly 
chosen spin. A value r = 10 will be found to be satisfactory. The new spin is then normalized 
to unit length. If the energy of the lattice is lowered by using this new spin, the new spin 
is adopted. Otherwise, the new spin is adopted with a probability exp[(E Q — E n )/T}. This 
move satisfies detailed balance, and so this Monte Carlo procedure will sample the Boltzmann 
distribution | 12| . A natural unit of equilibration time, the Monte Carlo step (MCS), is N 2 
iterations of this move on a. N x N lattice. Most runs last for 800000 MCS after an initial 
equilibration of 80000 MCS. For the case of A = 2, 3200000 MCS will be performed after 
an equilibration time of 320000 MCS. The properties of the lattice will be sampled every 50 
MCS. 

The number of disclinations can be counted by looking at all (N — l) 2 plaquettes on the 
lattice and determining whether a disclination of strength +1/2, —1/2, or is present at 
each plaquette. This determination is made by first defining r^-114 to be a counterclockwise 
ordering of the spins around the plaquette. Double-headed spins are converted into single- 
headed spins by = n 2 x — n 2 y and m iy = 2rii X n iy . The following angle is defined for the 
plaquette: 

t = 021 + 032 + 043 + 014 , (81) 

where 0^- = Q\ — 9j is constrained to be in the range (— 7r, tt), and 0j is the angle associated 
with spin rrij. The disclination density at this plaquette is defined by 

s = l- < 82 > 
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The disclination asymmetry is equal to the sum of this density over the entire lattice. By 
the analog of Green's theorem, it can be written as a sum of 8^ over the boundary of the 
lattice. 

The correlation length and the correlation time will be measured during a Monte Carlo 
run to monitor convergence. The position correlation function is defined as 

^(r)=iV- 2 5:{[n(x)-n(x-r)] 2 -l/2} , (83) 

X 

where the r and x are integer vectors on the N x N square lattice. Free boundary conditions 
imply that there will be significant effects of the boundary in this correlation function. The 
correlation length is defined in terms of the position correlation function by 

f = J dvrg(r)/ J dvg(r) . (84) 

Here |r| ranges from to 2 1 ^ 2 N. Fast Fourier transforms will be used to compute this 
quantity in 0(N 2 \n 2 N) time fL3 |. The time correlation function is defined as 

f(t)=J2 m (t')-m(t'-t) , (85) 
v 

where the m(t) = J2 X m ( x > t) is the vector sum of the single-headed spins at MCS t. A 
correlation time could be defined by analogy with Eq. fl8~4]), but noise in f(t) for large t 
causes this approach to be unsatisfactory. Instead, the correlation time, r, is definedby the 
smallest value of r for which /(t)//(0) < 1/e. 

The interesting observable is the asymmetry between the number of +1/2 and —1/2 
disclinations. This asymmetry will be small, and it will be important to quantify the statis- 
tical error in this observable. The difference Sn = ri+i/2 — n-1/2 must scale as N/£ by Green's 
theorem. The variance of this observable should scale with the observable and inversely with 
the number of independent samplings: 

-N/e 112 



(86) 



T/r 

where T is the total number of MCS. This variance is independent of A for small A, and so 
the constant c can be determined from the case where A = 0. In this limit, the values of r 
and £ can also be determined from the case A = 0. 

21 



Figure 5 shows the correlation length as a function of temperature for various lattice 
sizes. The correlation length grows at lower temperatures. It saturates at a fraction of 
the lattice size for even lower values of the temperature. The isotropic-ordered transition 
temperature can be identified by the inflection point of this curve. Finite-size scaling can 
be used to extrapolate the inflection point as a function of 1/N to iV — > oo. The result is 
approximately Jr/T = 4. Figure 6 shows the correlation time, in units of 50 MCS, as a 
function of temperature for various lattice sizes. This time also grows for lower temperatures 
and larger lattices. In fact, near the critical point, the expected scaling r ~ c£ 2 is observed. 

Figure 7 shows the observed number of disclinations as a function of temperature for 
various lattice sizes. As expected, the number of disclinations decreases with decreasing 
temperature. These curves should converge to a universal curve in the limit N — > oo. This 
curve is the total number of disclinations, both bound and unbound, so it does not go to 
zero at the critical point. The curves for different N, however, do intersect at a unique value 
of J/T. This value can be extrapolated as a function of 1/N to iV — > oo to determine the 
true critical point. The result is Jr/T = 3.4. This observable appears to produce a more 
reliable critical point than does the correlation length. 

Not shown are the correlation length and disclination densities for non-zero A. The 
dominant dependence on A was through A 2 . Non-zero values of A increased the density of 
disclinations and decreased the correlation length and time. Furthermore, analysis showed 
that the observed disclination density scaled like p ~ c£~ 2 , at least for large T, where the 
disclinations were unpaired. These qualitative effects are consistent with Eqs. fl46|) and ( p5| ) . 

The limitation of a finite-size lattice prevents the true scaling of the disclination asymme- 
try near the isotropic-ordered transition from being observed. One can, however, observe the 
asymmetry for these finite-size systems. Figure 8 shows the asymmetry function p+1/2 — P-1/2 
for A = 2 for various lattices sizes. Figure 9 shows the analogous results for A = 3. This 
observable can be extrapolated as a function of 1/N to N — > 00 to estimate the behavior for 
infinite systems. It is clear that the extrapolated result is zero, within statistical error. In 
fact, one can see that the disclination number asymmetry scales only with the circumference 
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of the system. 



VII. DISCUSSION 

We see that there is a statistical symmetry in this model. This symmetry enforces 
P+i/2 + P-i/2 — in the limit of a large system size, even in the presence of non-zero A. This 
symmetry arises because the disclination number asymmetry can be written as an integral 
of bounded terms over the periphery of the system, by Green's theorem. This relation, in 
turn, means that the disclination number asymmetry can scale at most with the linear size 
of the system. Indeed, this scaling was observed in the Monte Carlo calculations. With 
a definition of disclinations not susceptible to Green's theorem, this statistical symmetry 
would not be present. In this case, one would generally expect the disclination asymmetry 
to scale with N 2 /£ 2 instead of N/£. 

There are many parallels between the problem of a nematic liquid crystal with unequal 
Frank constants and a hexatic membrane fll^ , [T5|j . Disclinations mediate a melting transition 
in both cases, and the disclination free energy is logarithmic in the correlation length in both 
cases. For the nematic, unequal Frank constants cause the free energies of +1/2 and —1/2 
disclinations to differ. They only differ by a core energy, however, due to the renormalization 
of A. For the membrane, local buckling causes the the free energies of five- and sevenfold 
disclinations to differ. They, too, only differ by a core energy, due to renormalization of 
the membrane rigidity. As with nematic liquid crystals, there is a topological theorem 
that relates the number of five- and sevenfold disclinations to an integral over a boundary 
fT6|| . The natural way for this to happen is n.5 — 717 ~ cR/^a, where £a is now the hexatic 
correlation length. The prefactor again depends on the difference between exponentials 
of core energies. We conclude, then, that within the liquid phase of a membrane with 
hexatic symmetry, differing core energies simply renormalize the line tension and Gaussian 
rigidity. Near the liquid to hexatic transition, however, constraints imposed by a non- 
Euclidean membrane geometry frustrate the formation of hexatic order. The system can 
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relieve this frustration by flattening the membrane. One might expect, for example, that 
the preferred size of a vesicle undergoing a hexatic to liquid transition should scale like the 
hexatic correlation length. This conjecture is a worthy subject of future calculations. 
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FIGURE CAPTIONS 

Figure 1. 6 versus (p for a s = +1/2 disclination. 
Figure 2. 9 versus <fi for a s = — 1/2 disclination. 

Figure 3. E/J versus A for s = +1/2 (dashed line) and s = —1/2 (solid line). 
Figure 4. Definition of spins used in Hamiltonian ([3C|). 

Figure 5. The correlation length as a function of temperature for the case J = 3 and A = 
and N = 4 (dashed line), 8 (dot-dashed line), 16 (long-dashed line), and 32 (solid line). 

Figure 6. The correlation time as a function of temperature for the case J = 3 and A = 
and N = 4 (dashed line), 8 (dot-dashed line), 16 (long-dashed line), and 32 (solid line). 

Figure 7. The disclination density as a function of temperature for the case J = 3 and 
A = and N = A (dashed line), 8 (dot-dashed line), 16 (long-dashed line), and 32 (solid 
line). 

Figure 8. The disclination asymmetry as a function of temperature for the case J = 3 and 
A = 2 and N = 4 (dashed line), 8 (dot-dashed line), 16 (long-dashed line), and 32 (solid 
line). The error bars are ± one standard deviation. 

Figure 9. The disclination asymmetry as a function of temperature for the case J = 3 and 
A = 3 and N = 4 (dashed line), 8 (dot-dashed line), 16 (long-dashed line), and 32 (solid 
line). The error bars are ± one standard deviation. 
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